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Abstract: Exploiting the azimuthal angle dependence of the density matrices we 



construct observables that directly measure the spin of a heavy unstable particle. 
A novelty of the approach is that the analysis of the azimuthal angle dependence 
in a frame other than the usual helicity frame offers an independent cross-check 
on the extraction of the spin. Moreover, in some instances when the transverse 
polarisation tensor of highest rank is vanishing, for an accidental or dynamical reason, 
the standard azimuthal asymmetries vanish and would lead to a measurement with 
a wrong spin assignment. In a frame such as the one we construct, the correct spin 
assignment would however still be possible. The method gives direct information 
about the spin of the particle under consideration and the same event sample can 
be used to identify the spins of each particle in a decay chain. A drawback of the 
method is that it is instrumental only when the momenta of the test particle can be 
reconstructed. However we hope that it might still be of use in situations with only 
partial reconstruction. We also derive the conditions on the production and decay 
mechanisms for the spins, and hence the polarisations, to be measured at a collider 
experiment. As an example for the use of the method we consider the simultaneous 
reconstruction, at the partonic level, of the spin of both the top and the W in top 
pair production in e + e~ in the semi-leptonic channel. 
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1. Introduction 



The Standard Model (SM) of particle physics has been successful in explaining all 
the collider observables till date with a high degree of precision. This is remarkable 
considering that the particle content of the model is not complete since it requires 
a scalar spin-less particle, the Higgs. In the SM formulation, this particle is an 
integral ingredient of the mechanism of electro-weak (EW) symmetry breaking. This 
mechanism is still not well understood. For example, the fact that in the SM no 
symmetry protects the mass of the spin-less Higgs poses the hierarchy problem. 
Any solution to these issues brings in new particles and interactions at TeV scale 
with varying spin and gauge quantum number assignment. In a collider experiment, 
where almost all these particles are expected to be produced and decay to the light 
SM particles, the gauge quantum numbers, in principle, can be re-constructed by 
adding-up the gauge quantum numbers of all the observed light SM particles. Spin, 
on the other hand, shows up only in the distributions in various kinematic variables 
in the production and the decay sub-processes. Since the knowledge of the spin, 
along with the gauge quantum numbers, can enable us to distinguish amongst various 
candidate theories of physics beyond the SM (BSM) there has been growing interest in 
this subject recently in the context of the upcoming Large Hadron Collider (LHC) [[IJ 

i, 1 1, i, i, 0, i, i, gg, [n], m, m, m, m, m, mi, m m and ais ° in the context of 



proposed International Linear Collider (ILC) f2(| P, |2|, |j| H, gjj. 

Most of the BSM models have new particles that are partners of the SM particles 
based on their gauge quantum numbers assignments. However, the new particles may 
differ in the spin assignments. In models with supersymmetry (SUSY), the spin of 
the SM partner differs by | owing to the fermionic nature of the SUSY generators. 
There are however many other models such as UED where the spin of the partner is 
same as in the SM. In both kind of models a Z2 symmetry can be left over, leading 
to a heavy stable particle in the spectrum which can be the dark matter candidate. 
In SUSY models, the lightest neutralino, singlino, gravitino or axino can be stable, 
while in the models with universal extra dimensions (UED) the first Kaluza-Klein 
(KK) excitation of photon can be stable and are the dark matter candidate. These 
dark matter candidates can not be detected directly in collider experiments. Thus 
if these particles appear at the end of a decay chain, the re-construction of spin can 
be non-trivial, specially at a hadron collider like LHC. Moreover it would be also 
important to infer the spin, and other properties, of these dark matter candidates 
since these properties are important for the indirect detection of dark matter in 
astrophysical experiments. 

The spin of a (new) particle determines the Lorentz structure of its couplings 
with the other SM fermions and bosons. This, in a way, fixes its dominant production 
and decay mechanisms. In many cases a careful study of the energy dependence of the 
cross section around threshold can distinguish between spin-0 and spin-| particles 
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for example. Other methods to determine spin involve decay particles correlators. At 
the heart of these more direct methods is the decay helicity amplitude. For example, 
the helicity amplitude of a particle with spin-s and helicity A with — s < A < s 
decaying into two particles of spins si and s 2 , with helicity Zi j2 respectively, can be 
written as |E6| 



= \j 2 -^DZ{<j>A-4>)M s luh = ^^^e^dU0)Ml h , I = h-h. 

(1.1) 

Here M.f j is the reduced matrix element. This has been written most conveniently 
in the rest frame of the decaying particle. In fact the helicity here is the projection 
of the spin on the quantisation axis. The polar angle 9 is measured w.r.t. this quan- 
tisation axis and the azimuthal angle is measured around the same quantisation 
axis with freedom to chose the = plane. In most of the examples, it is useful 
to chose the production plane of the decaying particle as the = plane. Boosting 
along the quantisation axis will leave the value of the helicity unchanged. The an- 
gular distribution in these angles encodes the spin information through the rotation 
matrix D which factorises into an overall phase factor e l( - x ~ 1 ^ carrying the azimuthal 
angle dependence and the the d s xl function carrying the polar angle dependence. 
The latter can be expressed as [|7]] 



r S fff , _ V^W [(8 + X)\ (3 - X)\ (3 + l)\ (S-1)\} 1/2 

w;-2^ tl ( s + X -t)\ ( s -i-t)\ (t + i-\)\ 

X { C " S 2) { Sm 2) (L2) 

with — s < A, I < s. The sum is taken over all values of t which lead to non negative 
factorials. The differential rates have therefore polynomial dependence on cos 9 up to 
degree 2s and the azimuthal modulation coming from the off-diagonal elements of the 
density matrix ranges up to cos(2s0). One can construct observables to extract the 
degree of cos# and/or cos0 distribution. If the highest mode for, say, the azimuthal 
dependence cos(2s0) can be extracted this would be an unambiguous measure of the 
spin, s of the particle. 

Other methods have been used or advocated to determine spin. 

1. Exploiting the behaviour of the total cross-section at threshold for pair produc- 
tion [15|, ^2| or the threshold behaviour in the off-shell decay of the particle [0], 

2. distribution [|], |7|, [IT], [13], |22| in the production (polar) angle relying on a known 



production mechanism, 

3. comparing different spin assignments to intermediate particles in a process for 
a given collider signature % §, 0, §, 
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4. comparing SUSY vs UED for a given collider signature || |4], [| 11, 12, 14, 16 



21, |3, H, 25 



5. extracting the (cos#) 2s polar angle dependence |2|, [3l, ffj |8|, |14|, |20| , |2~2"1 , [2~8| , [29 



or 



cos 2s0 azimuthal angle dependence |23], |4|, Q of the decay distributions. 

The first four methods are indirect ways to assess the spin information subject to 
some assumptions and can only support or falsify a hypothesis. For example, the 
threshold behaviour depends not only on the spin but the parity of the particle as 
well U and for a particle of given spin it could be used to determine its parity []30| 



Further, it has been shown |22| that for pair production in an e + e~ annihilation, the 
threshold behaviour alone can not determine the spin of the particle. With f3 the 
velocity of the produced particle in the laboratory, at threshold the cross section for a 
scalar scales as /3 3 behaviour while for a spin-1/2 it goes like (3, except for Majorana 
fermions which can have /3 3 behaviour. Note that these (3 n characteristics do not 
take into account Sommerfeld / Coulomb f31|1 type corrections. The spin-1 particle 
can also have threshold behaviour and production angle dependence same as that of 
fermions with only difference coming from the distributions of their decay products. 
Thus, threshold behaviour and production angle distributions can at best be used 
only to confirm the spin assignment not to determine it. In the second method, one 
usually assumes a production topology, like for example s-channel pair production 
through a gauge boson. In this method the production angle dependence will depend 
upon the spin of the test particle. But still this dependence is not unique and can 
be obtained for higher spin test particles. 

The third and the fourth method uses numerical values of correlators or differ- 
ences in the distributions, which can be modified by the changes in the couplings or 
the presence of addition particles in virtual exchange etc. Thus, one can not use this 
method without having re-constructed the spectrum of the theory experimentally. 
The last method, which uses decay correlators, gives either the spin of the particle 
or the absolute lower limit on its spin. We note that the moments of the polar angle 



distribution discussed in Ref. [28] gives an upper limit on the spin of the particle. 



1.1 Spin through the polar angle 

Earlier studies of spin measurements used the average values of cos 9 or angular asym- 
metry, with appropriately defined polar angle 9, in the process of 2-body decay [^8 



or cascade decay f29fl . The numerical values of the angular asymmetries or the mo- 



ments of angular distribution gave estimates of the spin of the decaying particles in a 
model independent way. Most of the recent spin studies using decay kinematics focus 
on a decay chain that can be realised in SUSY or UED models. All the intermediate 
particles in the decay chains are assumed on-shell such that there is no distortions 
coming from the shape of the off-shell propagator and that it can be decomposed as 
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Figure 1: The typical decay chain studied for the spin analysis of particle A via its decay 
into observed particles B and D and a missing particle E. 



a series of two body decays, simplifying the calculations. For example, we look at a 
3-body decay chain of a particle A, shown in Fig.|I]. We look in the rest frame of the 
intermediate particle C whose spin is to be determined. Using crossing symmetry 



we write the matrix element for the s-channel process AB — > C — > DE as 



(9 BD ,(f>) = (2s + l) 



'BD ) 



;i.3) 



where, A, = X B — A^ and A/ = Ad — As. The rotation matrix d^. Xf (9) for spin s 
particle is 2s degree polynomial in cos(#/2) and sin(0/2), Eq. ( |1.2|) , which on squaring 
transforms to a 2s degree polynomial in cos 8. This leads to a 2s degree polynomial 
form of the angular distribution as 



dT(A -> BDE) 



Qo + Qi cos9 BD 



Q 



2.H 



COS 



2.s 



BD- 



;i-4) 



dcos9 BD 

Thus, we see that the degree of these polynomials is a consequence of the representa- 
tion of the particle under Lorentz or rotation group, in other words, the spin of the 
particle, provided C is produced on-shell, i.e. (p^ + Pe) 



Pc 



m 



c 



constant. 



One can also describe the decay in powers of some invariants, the highest power 



giving a measure of spin. Indeed, with m BD 
dT(A -> BDE) 



{Pb + Pd) 2 we can write 



dm\ D 



P + Pi m BD + ... + P 2s (m BD ) 



2.s 



1.5) 



obtained from Eq. (]1 .4]) through dm BD = 2E B E D f3 B f3 B d cos 9 BD by using a trans- 
formation of variables. Note that we could have P2 S = and P2S-1 7^ for some 
kinematical or dynamical reason, in this case we would set the lower limit on the spin 
to be s — |. We note that the above method involves two decay products of particle 
A while it measures the spin of the intermediate particle C and not the spin of the 
mother particle A. To directly measure the spin of A, we need to use the polar angle 
of every decay products w.r.t. the quantisation axis of A. The distribution w.r.t this 



decay angle looks identical to Eq.(O) with s being the spin of A. 
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1.2 Spin through azimuthal angle 

Another method of direct spin re-construction is to use the azimuthal angle distri- 
bution of the decay product about the quantisation axis of the decaying particle A. 
This is the main thrust of the present work. Using the form of the rotation matri- 
ces it can be shown, see later, that the azimuthal distribution appearing from the 
interference of different helicity states, has the general form 

dT 2s 2s 

JI = a o + ^2 a i cos( ^) + ^2 b i sin (i^)> ( L6 ) 

^ 3=1 3=1 

with a,j being the CP even contributions while the b 3 - being CP odd contributions. 
A statistically significant non-zero value of a 2s /ao or &2s/flo proves the particle spin 
to be s. The coefficients, a d - and bj, depend on the dynamics of production and 
decay processes and we will see that they are proportional to the degree of quantum 
interference of different helicity states of the particle A, or in other words, to the off- 
diagonal elements of production and decay density matrices. This distribution (but 



for the CP odd part) has been proposed in [£3| and used in [g4j to measure the spin of 
W and Z bosons at LEP-II and Tevatron, respectively. The azimuthal distribution in 
the laboratory frame is not simple sin or cos, however it is sensitive to the polarisation 
of the decaying particle as shown in Ref. []32]j . In this paper, we study the azimuthal 
distribution, Eq. (|1.6| ), in a model independent way to determine the constants ajs and 
bjS in terms of production and decay mechanism and construct collider observables 
to possibly measure these constants. We construct the observables in two different 
frames of reference and compare their merits. 

This paper is organised as follows. In section 2 we give the angular distribution 
of decay products for a general process of production and decay with emphasis on 
the case of spin— ^ and spin-1 particles. We describe the azimuthal distribution in 
terms of observables (asymmetries) to be used at colliders or event-generators in 
section 3. A numerical example of top quark decay chain is given in section 4 for the 
two different reference frames. We conclude in section 5. Additional expressions are 
given in the appendices. 



2. Density matrices, polarisation and azimuthal distributions 



To assess the spin of an unstable particle A, we look at a general n-body production 
process B\B 2 — > A A\ ... A n _i followed by the decay of A as A — > BC, for example. 
The other particles A^s produced in association with A can be either stable or decay 
inclusively. The differential rate for such a process is given by (see for example ||32||), 



A, A' 



(2nY 
21 



p(A, A')5 4 (^k Bl +k B , 2 -p A - Pi^j ^ 



d 3 PA 
2E A (2n) 



n-l 

n 



d?pi 
2Ei(2ir) 
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F~ o r <A A )<* (Pa - Pb - Pc) 



2m A v ' ' ^ A *° ^'2E b {2tjY 2E c (2tt) 3 _ 



:2.v, 



after using the narrow-width approximation for the unstable particle A, thereby 
factoring out the production part from the decay. Here we have I 2 = [m BlB — 
(m Bl + m B2 ) 2 }[m 2 BiB2 - (m Bl - m B2 ) 2 }, m| lBa = (k Bl + k B , 2 ) 2 , T A is the total decay 
width of A, rriA is the mass of A and Ta « rn A - The production and decay density 
matrices for A are denoted by p(X, A') and T'(A, A'), respectively. The terms in square 
brackets in Eq. (|2.1|) are Lorentz invariant combinations. The phase space integration 
can be performed in any frame of reference without loss of generality. 

Since we are interested in the decay distribution of A, we perform the phase space 
integrations in the rest frame A. We integrate the first square bracket in Eq. ( |2 . 1| ) 
and denote it as 

( 27r ) 4 f\ M\sif, , f \^ \\ ^?PA d 3 Pi 



^,A0=/^p(A,A0^(ife ft + ife*-px-(Eft))2^fe? II 



2E t (2n) 3 ■ 
(2.2) 

We note that the total integrated production cross-section, without cuts, of the 
process is given by the sum of diagonal terms a a = Tr <r(A, A'), while the off- 
diagonal terms of cx(A, A') denote the production rates for transverse/tensor polar- 
isation states or, in other words, for the quantum interference states. Further, we 
rewrite <r(A, A') = a a -Pa(A, A'), where Pa(A, A') is the polarisation density matrix for 
A in the corresponding production process. Similarly, we can partially integrate the 
second term in Eq. ([2.1|) and write it as 



d 3 p B d 3 p c 



2E B {2nf 2E c (2ir) 3 
B BC (2s+l 



An 



r A (x,\')dn B , (2.3) 



where B B c is the branching ratio for the decay A BC, s is spin of A, Ta(X, A') 
is the decay density matrix normalised to unit trace, dQ B is the solid angle measure 
for the decay product B. 1 Combining Eq. (|2.2j) and Eq.( |2.3| ) in Eq. ( |2.1| ) we get the 
decay angular distribution as 

-^f = ^E^( A ' A ') r ^' A '), (2-4) 

a A, A' 

where a = o a B B c, the total cross-section for production of A followed by its decay 
into BC state. The polarisation density matrix contains the dynamics of the pro- 
duction process and we will discuss its form for spin— ^ and spin-1 particle in the 



1 One can also consider 3-body or higher body decay of A in Eq.(2.1) and write Eq.(E^) by 



integrating all the phase space except f2 g . One example of this will be the top-quark decay j32[ 
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following sections. 



First we will discuss the general structure of the decay density matrix which can be 
studied independently of the production mechanism. The decay density matrix for 
a spin-s particle, expressed in terms of helicity amplitudes Eq. flLip , is given by 

r' s (A,A') = J2 M hiKt' 

h,h 

'2s + 1 



4:71 



e^- x '^J2 d xiWU(0) \M s hjl2 \ 



l L h 



2s + 1 
4tt 



\M s lull J 



where 



2s + 1 

4lX 



h 



\h\ < s i) | /i — / 1 < s 2 , \l\ < s. 



(2.5) 



(2.6) 



For a spin s particle there are 2s + 1 different af's that define the decay density 
matrix with 

Tr(r' s (A,A')) = E<" 
i 

Dividing the V s by its trace leaves us with 2s independent quantities involving af's to 
define the normalised decay density matrix of a spin-s particle. Now, the normalised 
decay density matrix can be written as 



r A (A,A') = e 



i(A-A') 



EX 



3 <(a - a/) *7a(A,A';0), (2.7) 



where 7a(A, A'; 6) = ja(X,X') is the reduced normalised decay density matrix with 
only 9 dependence left. It is important to keep in mind that the dependence is 
an overall phase and we see clearly that the differential cross section will have a 
more transparent dependence on the azimuthal angle than the polar angle. Using 
the relation in Eq. (|2.7|) we can re- write Eq.( [2.4|) as 

2s + 1 



1 da 



Pa{\ A) 7 a(A,A) 



a dQ,B 47r 

+ J2 &[*U(A,A')] 7a(A,A') cos((A - A')</>) 



A^A' 



%IPa(\\')) 7a(A,A') sin((A-A')0) 

A^A' 



(2- 
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which is similar to Eq.(L6) after integrating out cos (9. Thus we have a simple looking 
distribution of the decay product and also the coefficients of the different harmonics 
in the distribution. We emphasise again that the <fi dependence enters only through 
terms with A 7^ A', in other words the off-diagonal elements of the production and 
decay density matrices. When integrating over the full space without any cuts, the 
information contained in these terms will be lost. Another point to stress is that the 
form of the distribution Eq.( f2.8| ) remains same in any other frame as long as is 
measured around the momentum axis of the particle with some suitable reference for 
= 0. The measurement of the cosn0 (with n < 2s) modulation that stems from 
the part describing the decay depends on the size of the corresponding factor PaIa 
which is controlled by the interactions of particle A. The factors Pa describe the 
different polarisations with which the particle is produced and the factor 7^4 depends 
on the dynamics controlling the decay. One of the aims of this paper is to analyse 
how one can use this understanding to maximise these modulations and especially 
the modulation with cos 2s<p which is the most unambiguous measure of the spin-s of 
the decaying particle. For illustration and as a guide, in the following, we take a close 
look at the production and decay density matrices for spin-| and spin-1 particles 
to identify the conditions on the production and decay dynamics for the spin to be 
measured. 

2.1 Spin— I particle 

For the decay of spin-1 particle, — ► + | s 2 , ^2)5 the normalised decay 



density matrix, in the rest frame or rather the helicity rest frame [|33]], can be written 

as 



ri(A,A') 



l+a cos 9 a sin 9 icf> 

2 2 e 

a sin 9 p —i<f> 1—a cos 9 

2 C 2 



;2.9) 



using Eq.([2.5Q. Here a = (a^ — a -^ 2 /2)/( a i/2 + Q-i/z) anc ^ a \ are defined in terms of 
reduced matrix elements in Eq. ([2.6| ) and for spin-| particles given as 

= {£) E Kl-i /2 i 2 i^i < *i> ft - V2i < s 2 



1/2 
a -l/2 



(^ElKl+i/sl 2 I'll^L Hi + 1/2|<* 2 . (2-10) 



The explicit calculation of a for the spin-^ particle decaying into two body final 
state is given in the Appendices |C.1| and |C.2| for decays into a lighter spin-1/2 and 



either a scalar or spin-1 considering general effective operators. We have restricted 
ourselves to operators of dimension 4. It can be seen that a is zero for a pure vector 



9 



or pure axial- vector coupling in the case of decay to a spin-1. It can also be small 
depending on the masses of the daughter particles. 

The polarisation density matrix for a spin-1 particle can be parameterised as 



Pi (A, AO 



1+^3 Vl- ir l2 

rji + ir] 2 1 - i] 3 



(2.11) 



where rji is the transverse polarisation of A in the production plane, 772 is the trans- 
verse polarisation of A normal to the production plane and i] 3 is the average helicity 
or polarisation along the momentum of A or polarisation along the quantisation axis. 
Combining the expression of Ti(A, AO and Pi (A, AO in Eq. (|2.8|) we get the angular 
distribution of spin-| particle as [fffil 



1 da 1 1 ,. . . . . . n • 

— — = — 11 + ar] 3 cos u + am sin U cos + ar] 2 sin v sm . 

The cos 9 averaged azimuthal distribution is given by 



1 da 1 

0"i d(f) 



1 

2^ 



ar/iir 

1 H cos ( 

4 



ar] 2 n . 
H — sm ( 



(2.12) 



(2.13) 



Here we note that the cos0 or the sin0 modulation of the azimuthal distribution 
is proportional to the transverse polarisation of the spin-^ particles and also to 
the analysing power a. Thus, it is important that the production process yields 
a non-zero value of either 771 or ?y 2 - A non-zero r/ 2 indicates CP-violation or the 
presence final state interaction (absorptive parts). A non-zero r/i can be obtained 
either with parity violation, which is present in the electro-weak sector of the SM 
or with appropriate initial beam polarisations. Further, we also need to know the 



analysing power a of the particle. For spin-^ particle, it is given in Eqs. 



and 



( |C.6| ). We see that the decay vertex has to be at least partially chiral, i.e. parity 
violating for a 7^ 0. That is, we need effectively chiral production and at least 
partially chiral decay for the fermions for their spin to be measured. 
It is educative to realise that Eq.( |2.12| ) can be cast into 



1 da 1 
0i dfl 



B 



1 

47T 



1 + a 



Pb 



■n 



B\ 



with if = (771,772,773) 



(2.14) 



with 77 the polarisation vector. Performing a general rotation will leave pb-V un- 
changed. In the new frame, after rotation, we can define a new averaged azimuthal 
distribution as 

1 da 1 



1 

2^ 



1 H 9— cos ( 

4 



/ H — sm . 

4 



(2.15) 



If the rotation is done along the r\ 2 direction (normal to the production plane), then 
r] 2 = rj2 but 7][ will pick up a contribution from 773, the average helicity. If r] 3 ^> 771 the 



10 



azimuthal distribution in this new frame is more conducive to a spin measurement, 
in the sense of catching the cos dependence. 

It is important to observe that the picture we have described so far in terms of 
azimuthal dependence through cos and sin (or higher for higher spins) may be very 
much impacted if cuts are applied to the cross section. If the cuts are 0-dependent, 
the azimuthal distributions may no longer have the simple form of Eq. ( [2.1 3D but 
would carry "spurious" dependence that would prevent the spin reconstruction as 
advocated here. Indeed, we could have a much more complicated dependence of the 
form 

1 d<Ji 



0"! 



1 
2^ 



Feit) + ^G c (0) cos + F c (0) sin , 



(2.16) 



unless only independent cuts are applied as suggested in Ref. [E!J 



2.2 Spin— 1 particle 

For the decay of spin-1 particle, |1, /) — > h) + \s2, ^2), the normalised decay density 
matrix is given by 



ri(U') 



where, 



and 



l+<H-(l-3<5) cos 2 9+2a cos 9 sin 6>(q+(l-3<5) cos 9) i q 

4 2V2 e 

sin6>(a+(l-35) cosfl) g -i</> ^ _|_ ^ _ 3^ sin 2 9 

(l_ 3( 5)(iz£2!! 



;i - 35) 



(l-cos 2 6>) j20 



sin0(a— (1— 35) cos0) j0 
273 e 

_i20 sin 9{a- (1-35) cos 9) -iff, l+<5+(l-3<5) cos 2 6>-2acos B 



2V2 



a] — a 1 



a 



a \ + a l + a -i 



v 7 h 
7 h 

(|)ek,, 1+1 i 



a l + a + a -l 

\h\ < 8 X , \h - 1| < s 2 

< min(si, s 2 ) 
I /1 1 < Si, |/i + 1| < s 2 



(2.17) 
(2.18) 



(2.19) 



The explicit calculation of the analysing power parameter a (the vector part) and 
5 (the rank-2 tensor) for a spin-1 particle decaying in two body final state is given 
in the Appendices |C.3| , |C.4| and |C.5| . In particular 5 = for decays into massless 
fermions assuming dimension-4operators. For the decay W — > ff in the SM we 
have a = — 1 for massless / and /'. 
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The polarisation density matrix of spin-1 particle has two parts: the vector 
polarisation which we define here as p = {p x ,Py,Pz) and is identical to that for a 
spin-i ff and the tensor polarisation described through a symmetric traceless rank-2 
tensor T^, Tr T=0. The density matrix is parameterised as |3c 



•Pi (A, A' 



1 _i_ 'pz_ _i_ T zz Px ip-y I T X z iTy Z 

3^2^^ 2^ %/3 



2^ V3 



1 _ 27^ 

3 VE 



T X x Tyy-\-1iT X y Px~\~ipy T X z~\~iTy 



T<£ X Ty y HT X y 

Ve 

Px iPy Txz ^Tyz 

1 Pz I T zz 
o 



(2.20) 



VE 2^2 V3 3 2 1 v^S 

Again using Eq. (|2.8|) we can write the angular distribution for spin-1 particle as 



1 dcr 2 3 

<T 2 dVL B 87T 



I - (1-3(5) ^ + ap z cos# + ^(l-35) T zz cos 2 6 



+ yOL p x + 2y ^(l — 35) T xz cos9j sin6 cos( 
+ ^ct p y + 2y^(l — 35) Tj^cos^ sin 6* sin<j 



+ 



1 - 35) ( 



T 

J- .7' 



35 ) ) sin 2g cos(20) 



h — 35) T xy sin 2 fl sin(2</>) 



(2.21) 



The cos # averaged distribution is 
1 da 2 



0-1 



3 

47T 



2 ap^vr ap v TT 2 

- + — ^— COS0+ — ^— S1B0+ -(1 



xy , sin(20) 



35) 



( 



T -T 



V Vq 




"" 1 .-<»s(2o) 

(2.22) 



We note that the 4> and 2(f) modulation of the azimuthal distribution is proportional 
to the transverse polarisations, p x and p y , and the transverse components of the 
tensor polarisations, T xx — T yy and T xy . Again, in this frame all the cos modulations 
are CP-even and the sin modulations are CP-odd. To determine the spin we need 
T xx — T yy 7^ from the production part and 5^1/3 from the decay part in the CP- 
even production process. Since there is no symmetry that sets 5 to be one-third, the 
dynamics in the decay part is not constrained. In other words, the decay mechanism 
does not require any parity violation. 

As we have done for spin-1/2 it is instructive to rewrite Eq.( |2.21[) in terms of 
invariants under rotations. If one defines a rank-2 tensor out of the tensor product 
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of the unit vector describing the momentum of the decay product B } ¥b = Pb ® 
Pb/\pb\ 2 2 , using the fact that T is traceless Eq. Q2.21p writes in terms of (rotation) 
invariants as 



1 dao 1 



a 2 cLVLb An 



1+a lm p+(1 - 3l) n TrB 



(2.23) 



We can then rewrite Eq.( |2.23|) in another frame, in particular one where we make a 
rotation around the y axis, transverse to the production plane. This will not mix the 
CP-odd and CP-even tensors but may make some <ft asymmetries in the new frame 
larger. 

2.3 spin— | and spin— 2 

For spin-| and spin-2 particles we give the decay density matrix in Appendix |B|. 
Since we need the coefficient of the highest harmonics to be non-zero for the spin to be 
determined, we note that for spin-| we need parity violating interaction in the decay 
process, i.e., ot\ ^ and/or a 2 ^ whose combination defines the analysing power 
of highest rank. For spin-2 particles we need A 4 oc (a| — Aa\ + 6oq — 4a' 2 i l + a?_ 2 ) 7^ 0. 
A4 is the analysing power of rank-4, the highest rank for spin-2 to be non- vanishing 
(see Eq. QB.5Q ). This can be achieved without parity violating interactions in the 
decay process. We note that parity violating interactions are required in the decay 
of fermions for its spin to be measured along with its (transverse) polarisation being 
non- zero . For the bosons, on the other hand, we only need its transverse polarisation 
being non-zero either due to parity violation in the production process or due to 
polarisation of the initial beams. 



3. The azimuthal distribution at event-generators/colliders 



The azimuthal distributions Eqs.( [2.13| ), ( [2.22 ) etc. are given in the rest frame of the 



decaying particle. To be able to measure the spin we need to construct the above 
mentioned azimuthal angle in terms of quantities defined in the lab frame of a collider 
experiment. Before considering other frames let us first define some asymmetries in 
the rest frame. 

3.1 Asymmetries in the rest frame 

We start with re-writing the rest frame azimuthal distribution in terms of some 
simple asymmetries that we define below. Let us first define 



(3-1) 



''(PB ® p B )ij = PB i PB j • The scalar product is TP B = £V ■ T„P B ij = Tr TP B . 
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For s = 1/2 we define the following asymmetries and calculate them using Eq. (|2.13 ): 

/i(-7r/2, tt/2) - Ji(tt/2, 3tt/2) cwji 



Ii(0,27r) 

7l(0,7r) - Ji(7T, 27r) _ m/2 

/i(0,2tt) 2 ' 



(3.2) 



These asymmetries have been used in Ref. |32[ as a probe of the polarisation of the 
top-quark. Eq. ( [2.13 ) can be re- written in terms of these asymmetries as 



1 dai 1 

~ 2^ 



7rA\ , TlB\ 

1 H cos d> H sin < 



(3.3) 



Similarly for s = 1 we further define similar asymmetries and calculate them in terms 
of vector and tensor polarisations as follows, 

/ 2 (-tt/2, tt/2) - / 2 (tt/2, 3tt/2) 3ap x 



A " h(0,2n) 4 

B x = h(0, 7r) - /2(7r,27r) = 3ap y 
2 ~ / 2 (0,2tt) ~ 4 

J 2 (-7r/4, tt/4) - J 2 (7r/4, 3tt/4) + J 2 (37r/4, 5vr/4) - / 2 (5tt/4, 7ti/A) 



Al 



J 2 (0,27r) 



2 n_ 3( j) 

V \/6 



5 2 2 



/ 2 (0, tt/2) - / 2 (tt/2, tt) + J 2 (7r, 3vr/2) - / 2 (3tt/2, 2tt) 



I 2 (0,27r) 



z(l-3*)f\/|^ 



7T 




(3.4) 



Thus Eq.( p.22|) can be re- written as 
1 d<j 2 



a 2 



1 

2^ 



nAl 



COS0 + 



sin + 



cos(20) 



txBI 



sin(20) 



(3.5) 



We see that the <fi distribution in the rest frame of the decaying particle has very 
simple form in terms of the above mentioned asymmetries. It is clear that for higher 
spins we need to cut the 2tt in more and more parts. The important observation 
to make is that the coefficient of the cos0, A\ for spin-1/2 and A\ for spin-1 are 
determined exactly in the same way in terms of the asymmetries, this generalises 
also in the same to higher spin particles and similarly for other coefficients of cos j(f>. 

Next we write the asymmetries in terms of spin-momentum correlators. To this 
effect, we first define the following spin vectors in the helicity rest frame 



s x = (0,1,0,0), s y = (0,0,1,0), ^ = (0,0,0,1) 



(3.6) 
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Modulations 


Asymmetries 


Spin-momentum correlators Cj 


cos(0) 


A 1 


Sx-VB 


sin(</>) 


B 1 


S y .p B 


cos(20) 


A 2 


(Sx-PB) 2 - (Sy.p B ) 2 


sin(20) 


B 2 


{Sx-PB){Sy.p B ) 


cos(30) 


A 3 


(sx-Pb) 3 ~ 3(s x .p B )(s y .p B ) 2 


sin(30) 


B 3 


3(S x .p B ) 2 (Sy.p B ) - (Sy.p B ) 3 


cos (40) 


A 4 


{Sx-PBY - 6(S X .p B ) 2 (Sy.p B ) 2 + (Sy.p B ) 4 


sin(40) 


B A 


(S X -PB) 3 (Sy.p B ) - (S X -P B )(Sy.p B ) 3 



Table 1: The table of asymmetries A 3 and B J corresponding to the j(ft modulation of the 
azimuthal distribution and corresponding spin- momentum correlators Cj. Here s xy are the 
transverse spin directions of the decaying particle A, with s x being in the production plane 
and p B is the 4-momentum of the decay product B. The spin vectors s x , y are listed in 
Table |2| in different frames. For frame F we replace Si with §i in the above correlators. 



which are orthogonal to the 4-momenta of the particle A, p^ = {rriA, 0, 0, 0). These 
spin vectors satisfy the conditions p.Si = and Si.Sj = —Sy. The asymmetries for 
the spin-i case can then be written as 

A x = v(s x -Pb < 0) - a(s x .p B > 0) B i = (y{sy.p B < 0) - cr(s y .p B > 0) 
1 cr(s x .p B < 0) + a(s x .p B > 0) ' 1 a(s y .p B < 0) + cr(s y .p B > 0) ' 

The asymmetries for the spin-1 case can be written as 

! = a(s x .p B < 0) - o{s x .p B > 0) ! = <r(s y .p B < 0) - a(s y .p B > 0) 

a(s x .p B < 0) + a(s x .p B > 0) ' 2 a(s y .p B < 0) + a(s y .p B > 0) ' 



A\ 



v({Sx-Pb] 2 - [SyPB] 2 > 0) ~ a([s x .p B ] 2 - [SyPB] 2 < 0) 
<t([Sx-Pb} 2 - [SyPB? > 0) + a{[s x .p B } 2 - [SyPB? < 0) ' 



B 2 = <7([Sx-PB][Sy.PB] > 0) ~ (?{[Sx-PB\[SyPB\ < 0) , g ^ 

2 <K[s x -Pb][VPb] > 0) + <r([s a .p fl ][srP B ] < 0) 



It is clear from Eqs. (|3.7|) and Q3.8|) that the asymmetry corresponding to a given 
modulation of </> has identical expressions in terms of the spin-momentum correlator 
s x -Pb and s y .p B . What we mean is that by constructing specific functions with 
products of Si.p B one reconstructs the set of cosm0, sinm0, see Table [I]. Thus in a 
spin independent way we can write these asymmetries as 

A > or B > = (_!,, °fj > ^ - ggj < °> , (3.9) 
v ; a(Cj > 0) + tr(Cj < 0) ' v ' 
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Figure 2: The momentum configuration in the laboratory frame L of the colliding beams 
is depicted. The angle between the production plane (the xz-plane) and the decay plane 
(spanned by pb and pc) is the azimuthal angle 4> M (= </>|/), which is the ordinary azimuthal 
angle in the frame M. This azimuthal angle has been studied in [23| and pi]. 



where the correlators CjS are listed in Table [I] for different modulations. Further we 
note that these expressions of asymmetries have a simple interpretation in terms of 
the polarisation parameters as long as s* z , 3- vector, is parallel to the 3-momentum 
of the decaying particle in the frame of choice and Si are orthogonal to each other. 
This defines a helicity frame. The lab frame, achieved by a boost along z-axis and 
then rotation around y-axis, also satisfies the properties of being a helicity frame. 
We note that the orthogonality of Si is respected only if the boost is along one of 
Si directions. In some other frame where these properties are not valid one needs 
to re-write these s^s as linear combinations of orthogonal Si as we will see in the 
following sections. This will be necessary if the new frame is not reached through a 
boost made along the direction of motion of the particle. 

3.2 The rotated frame M 

The rotated frame is in fact the frame that is obtained from the rest frame by 
performing a pure Lorentz boost along the quantisation axis, the amount of boost is 
such that the energy of the particle whose spin we want to study is the same as the 
one measured in the laboratory frame. This is therefore a helicity frame in the sense 
that the quantisation axis has now been identified to lie along the momentum of 
the particle. The normalised 3-spin vectors si remain unchanged and therefore also 
the polarisation vectors (fj, p) and other tensor polarisations. Hence the azimuthal 
asymmetries are the same as in the rest frame and have a one-to-one correspondance 
to the polarisation tensors defined in the rest frame. The appellation rotated frame 
comes from how this frame is picture in the laboratory frame. In fact this can be 
viewed as a simple rotation. In the laboratory frame L, defined in Table 0, the 
production plane of the particle A defines the xz-plane and the plane containing the 
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decay products with momenta ps and pc defines the decay plane. These two planes 
intersect along the momentum pa of the decaying particle, see Fig. |2j. Thus the 
angle between these two planes is the azimuthal angle of the decay product around 
the axis of spin quantisation (the momentum p^) , i- e. the we have mentioned in 
Eqs.( [2.13"D and ( |2.22| ) and which corresponds to exactly the azimuthal angle defined 



in the rest frame. In terms of the variables defined in the laboratory frame L, it is 
defined as 



s M .pg s L .pg 
where, 0™ f = tan -1 -|r — rr = tan^ 1 — — r . (3.11) 

Sx-Pb s x-Pb 



From Eq. (|3~T0|) it appears that the reconstruction of this angle in the laboratory 



requires that one measures the momenta of all the decay products. This may be 
achieved even if C, say, is invisible provided one has enough control and constraints 
on the initial state so that the momentum of the decaying particle A is known, like 
for instance in e + e~ annihilation where the beam energy is fixed. 

Another way to get the azimuthal angle 0a/ it to re-construct the scattering 
angle 9\ in the lab frame and rotate the event about the y-axis by that angle to 
bring the momentum p& in the direction of the z-axis. The azimuthal angle of the 
decay product, 0|f is same as mentioned above in the lab frame. We dub this frame 
the rotated frame M (obtained by rotating the laboratory frame) and the momenta 
in this frame as compared to that in the laboratory frame is given in the Table g. 
Using the form of the momenta ps and Si in frames _R(rest frame) and M, one can 
see that = 0# = 0|f. Thus the distribution in angle 0, defined in Eq. (|3.10| ), is 
the same azimuthal distribution as in the rest frame with same amplitudes for the 
different harmonics. 



This azimuthal angle has been first studied in Ref. p3 to demonstrate the simple 



cos(j'0) modulations of the azimuthal distributions and has been used to examine the 
spin of Z and W bosons at Tevatron and LEP-II, respectively, in Ref. [p^j . Here we 
provide a theoretical understanding of the amplitude of these cos(j0) modulations in 
terms of transverse polarisations of the particle under consideration and its analysing 
power a etc. . In the event when the transverse polarisation is negligibly small or zero, 
this frame will not give any modulation in the azimuthal distribution as shown in 
section |] for top pair production in e + e~ . To address this potential issue we construct 
another frame which will give us an independent estimate on the modulations of the 
distribution and hence the spin of the particle. A hint on how to achieve this has 
been illustrated in section 2.1 and section 2.2 for the spin-1/2 and spin-1 when a 
simple rotation mixed the longitudinal polarisation and the transverse polarisation 
in the production plane, leaving the polarisation transverse to the production plane 
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Rest frame := R 



(0,1,0,0) 
(0,0,1,0) 
(0,0,0,1) 



p* = (m A , 0,0,0) 



Pb 



E R 

B 



( 1 

f3§ sin 9§ cos 0^ 
P§ sin 9§ sin 0^ 
\ /3fcos^ 



Rotated frame := M 



(0,1,0,0) 



sf = (0,0,1,0) 



.s-." =(/?i,o,o,i) 7 i 
1, 



p^ = Ei(l,0,0,/?i 



V 



/3|sin< cos0f 
P%mn9% sin0^ 
(3% cos 6% 



J 



Lab frame := L 



(0, cos #i,0,- sin 0jj) 



^ = (0,0, 1,0) 
4 = (/3i,sin^,0,cos^) 7 i 
V \ = E L A (l,[3 L A sin 0^,0,/^ cos 0£) 
/ 1 \ 

sin cos 05 
sin 6>g sin 0^ 
/3g cos 



Pb 



E L B 



Boosted frame := F 



F 

Pa 



F 

Pb 



(0,1,0,0) 
(0,0,1,0) 
(^,0,0,1)71 



^(1,0,0,/?!) 

/ 1 

/3fsin0£cos0g 
j3g sin 9g sin 0g 
Pb cos 6*^ 



not the result of 
boost L — > F, 
see text. 



V 



Table 2: Momentum pa, pb and the spin directions Si in various frames. The transforma- 
tion R — > M is a boost along z-axis A 2 (/3^), M — » L is a rotation R(6 A ) and &: p|J are 
obtained by a boost along x-axis A x (— /3^sin0^) from frame L. Note that the spin vectors 
sf in frame F are not related to sf through boost but constructed such that they represent 
the helicity basis. The expressions for sf = A x (—P A sm9 A )sf which are the result of the 
boost are given in Eq. |3.12 . The azimuthal angle of interest is = tan~ 1 (s y .pB / s x .pb) in 
each frame (with s, replaced by in frame F). 



unchanged. The next section will show how this can be achieved in general and how 
we can construct the spin-momentum correlators in this case. 

3.3 The boosted frame F 

The idea behind the boosted frame F is to induce a non zero azimuthal asymmetry 
even in the event that transverse tensor polarisations are very small or vanishing by 
making the longitudinal components, assuming it is non zero, contribute. We will 
show how this can be achieved especially how we can construct the correlators from 
combinations of variables measured in the laboratory frame. It should be added 
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ny 



x' 




Figure 3: The momentum configuration in the transversely boosted frame F is depicted. 
The azimuthal angle of ps is the azimuthal angle <f> F (= </>f )• This frame is obtained from 
frame L by boosting along negative x-axis such that the transverse momentum p A becomes 
zero. 



that both the rest frame and the M frame are helicity frames. In the new frame 
and in order to arrive at the mixing between the longitudinal and the transverse 
polarisations we need to perform a transformation that will move the longitudinal 
spin (quantisation axis) away from the momentum of the particle. Yet, we still need 
to reconstruct a helicity basis in order to construct the helicity density matrix. To 
achieve the misalignment, we observe that a Lorentz boost in a direction other than 
the direction of the particle momentum will mix the helicity states. In our case one 
way to achieve this is to carry a boost from the laboratory frame along the negative 
x-axis with velocity @ A 'sm9 A ', thus reaching the boosted frame F. The momentum 
p F A of A is then pointing along the z-axis, see Table 0. This looks as if we have 
slowed down the particle, however contrary to frame M where the momentum is also 
pointing in the z direction, we can check that none of the transformed spin vectors 
sf = Aj.(— sf has its three momentum lying on the momentum of the 
particle. In fact, it can be shown on general ground that if initially the spin axis is 
parallel to the particle momentum, in the new frame these two axes will move away 
by an angle u, the Wick angle fl33| , if the boost is not performed along the particle 



momentum. Indeed we verify that 

sf = A x (—f3 L sm9 A ) sf = cosu; sf — since* sf 

Oy b y b y 

sf = A x (— (3 A sin 9 A ) sf = sin a; sf + cos lu sf with 

cos 9h . sin 9h /0 ^ nX 

cosu; = — , sinw = , (6.12) 

v/l-(^sin^)2 tV1-(#l sin 6^)2 

if y z are the helicity basis in the new frame F and are given explicitly in Table ^| 
Since the spin vectors sf z in the laboratory frame are not parallel to the x-axis 
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(except for sin 6^ = or 1), a boost along the x-axis modifies the orthogonality of 
the spatial component of Sj in the boosted frame. Also the spatial components of 
the boosted s z , i.e. A(— P^sin8^)s^, is not parallel to the spatial component of p F , 
owing to the Wick rotation of the spin basis. Thus the definition of the longitudinal 
or transverse polarisations in the frame F, which is not a helicity frame, is different 
from that in the helicity lab frame L. Since the asymmetries and in the frame 
F are defined with respect to sf, we can in principle have A 7 non-zero even in the 
absence of any transverse polarisation in frame L. The helicity basis s F y z is to be 
used to construct the spin-momentum correlators in frame F. 

This rotation of the spin basis vectors leads to the transformation of the density 
matrix and the various polarisation parameters. In general, for a rotation defined 



through the, Euler, angle 8, <f) the density matrix transforms as [33[] 

f/ = D$,9,-$) P M D\4>,e,-4>) . 

In our case we have 8 = u and <p = 0(boost in the x direction), which leads to 

p F (\,\') = d s M (u) d s xll ,(8 u ) p L (l,l') . (3.13) 

Thus the density matrix p F does not receive any additional phase and the azimuthal 
sin(n0^) dependence remains unaltered. The polarisation parameters, for the spin-1 
for example as concerns the vector and the tensor polarization, transform as 

pf = pf , T F = R lk {u) R jk (u) 2j (3.14) 

as mentioned earlier in Sec. [2.1| and [2.2[ Here, k, I = {x, y, z}, Rij is the matrix, 
which for the boost we have performed, corresponds to a rotation about the y-axis 
in the cartesian coordinate. The superscript F or L stands for the quantities defined 
in frame F or frame L respectively. The asymmetries and B-i in frame F have 
exactly the same expression as Eqs ([3.2|) and ( p.4[ ) with pi and Tjj replaced by the 
ones defined in Eq. ( |3.14| ). Thus for a spin-1 particle we have A 1 oc p F = R X j(u>)pj 
and not simply related to the transverse polarisation p L as we have in the (helicity) 
frames M or L. This shows again that one can have a non-zero A 1 even in the 
absence of transverse polarisation in the frames M or L. 

The azimuthal angle in this boosted frame is denoted by 0^ and shown in Fig.[3]. 
Since 0^ is the azimuthal angle around the new momentum p F Al it will have a simple 
cos(j0) and sin(j0) modulations in the distribution up to j = 2s. The asymmetries 
A 7 and B-i in this frame are defined w.r.t. the spin directions s F given in the Table 0. 
0^ is expressed as 

, ,f Sx-Pb cos^ s L x .p L B + sinw s L z .p L B 
cot <!>b = ^F-F = TlZZ ( 3 - 15 ) 

b y -Pb b y -Pb 

= coscu cot0^ + sincj %^f. (3.16) 

Sy-P B 
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We see that 0# is related to in a non-trivial way and thus the corresponding 
modulation need not have vanishing amplitudes even when this is the case in terms 
of distribution. Note that the asymmetries can be zero in both frames if either 
we have Pa(X, X') oc 5\ t \/, i.e. when the particle is completely unpolarised, or when 
the particle is spin-0. In all other cases, the two frames will lead to different values 
of the azimuthal asymmetries. Thus, we need to use both frames to confirm the spin 
of the particle. 

3.4 Note on event reconstruction 

The asymmetries A 1 and and the azimuthal angles <p M and (j) F require complete 
reconstruction of the test particle's momentum in order to construct the correspond- 
ing spin vectors sf and/or s^ F . The possibility of reconstruction depend both on 
the kind of collider and the number of missing particles in the process. For example, 
reconstruction of spin vectors is possible at colliders with fixed center of mass energy, 
like the ILC, for some selected processes where the number of missing particles is 2 
or less. At hadronic colliders, having variable centre of mass energy at partonic level, 
such reconstructions can be achieved for processes with one or no missing particles. 
Most of the new physics models with a dark matter candidate have two missing 
particles in the production process of new particles at LHC This makes the desired 
re-construction as outlined here unfortunately impossible at LHC in such processes. 
It is worth further investigating how this method could be combined with other 
methods or improved. To illustrate the method we therefore turn to an application 
for a collider such as the ILC. 

4. Application to e + e~ — > tt — > bW + bW~ — > bl + v bjj 

In this section we study top-quark pair production in e + e" in the semi-leptonic chan- 
nel as a test bed for the spin measurement based on the exploitation of the azimuthal 
asymmetries in different frames outlined previously. We chose this particular process 
because it represents a decay chain where the intermediate W boson is on-shell. The 
charged lepton, / + = e + ,/z + , in the leptonic decay of the top (and the W + ) will 
play the role of our particle B in the previous section and used to construct the 
spin-momentum correlators. In this example all the momenta can be reconstructed 
and therefore the methods we have outlined can be applied readily. We do not take 
beamsstrahlung into account nor do we consider the issue of backgrounds that might 
force us to introduce cuts, which we want to avoid. However, we consider the effect of 
beam polarisation. The polarisations of the initial electron and position beams can 
be used to tune the polarisation of the produced heavy particles that can drastically 
affect the polarisation. We work at y/s = 500GeV where the total cross section, 
including branching ratios, is 81fb for unpolarised e + , e~, This corresponds to a total 
number of 40500 events with a typical luminosity of 500 fb _1 . For each fit we make, 



21 



we will indicate what the minimum number of events is required for a 3a discovery 
of a particular cosj4> modulation that is a measure of the spin of the particle, we 
will see that this programme could be successfully carried at a linear collider with 
500 fb _1 for this process. 

For event generation we use the partonic level event generator Pandora-2 . 3 [[34] 
and generate 2 10 6 events for different initial state polarisations. Event by event we 
need to calculate s^f.p^, sff .p B (frame M) and s^.p B , Sy.p B (frame F) which can 
be expressed in terms of energies and angles measured in the lab frame, 



. p M = s L x .p L B = -£^(cos6^sin6)£cos0! - sin 6 L A cos 6g) 

T?L( oh ■ oh _ ■ oL „» aL\ 



~F F _ E b(I3a sin °A ~ sin °B COS 0|) 

SX ' PB ~ v/l-(^sin^)2 : 

»?-P% = = Sy-PB = ^-Pb = ~E L B sin 9 L B sin ^ . (4.1) 

We then calculate, for A = t,W all the 8 asymmetries corresponding to the correla- 
tors in Table [l], therefore testing whether a value for the top spin as high as s = 2 is 
possible. The azimuthal angles in frames M and F can be constructed using Eq . (|4 . 1| ) 
along with Eqs. ( |3. 1 If) and (|3.15|) for generating the distributions. The reconstructed 



azimuthal distributions are then fitted with a general function 

n 

F n {(j>) = a + J2 l a J cos(j0) + bj sin(j»] (4.2) 

3=1 

with n — 4. With n = 4, the only bias is that the particle has spin s < 2. We then 
compare the best fit coefficients with the asymmetries calculated. Since we work with 
the SM production and decay mechanisms for the t-quark, there is no CP violation 
in this process. The fitting procedure returns bj m and E>i pa in both frames M 
and F for all the initial state polarisations. This constitutes therefore a consistency 
check and confirms the absence of CP violation. In the following sections we will 
only talk about the CP even contributions coming from various cos(j'0) modulations 
and ignore the discussion on sin(j0) modulations as they are zero. 

4.1 Spin— | case: t— quark 

The top pair production at an e + e~ collider proceeds through a photon and a Z- 
boson exchange in the s-channel. We will study the effect of the initial polarisation of 
the electron P e - and positron P e + . The partial chiral nature of the Z coupling leads 
to a finite top polarisation even for unpolarised initial state electron and positron 
beams. For t-quark decaying into a lepton through a W, the analysing power of the 
top is a = 1. 

We start our analysis with unpolarised beams and the polarisation of top for this 
case is given as (P eT is the polarisation of e T ): 

(P e -,P e +) = (0.00, 0.00) : 77x = +0.222, r/ 2 = 0.000, r? 3 = -0.127 . 
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Quantities 


Frame M 


Frame F 


Reference 


(+0.00, +0.00) 


A 1 


+0.111 


-0.035 








— 


ai/a 


+0.175 


-0.055 






A 1 


-0.253 


+0.149 


Fig-! 


(+0.80,-0.60) 


ai/ao 




-0.397 


+0.234 


(+0.792, +0.60) 


A 1 


» 


+0.021 


Fig.| 






ai/a 


» 


+0.033 



Table 3: Values for the fitted asymmetry A 1 and the fit parameter ai/ao, see Eq. 4.2 , for 
the lepton distribution from i-decay for different initial state polarisations P e - , P e + for the 
electron and the positron in frames M and F. We have the relation ai/ao = nA 1 /2 which 
is observed numerically within tolerance (+10 -3 ). The other aj/ao,j ^ 1 are zero within 
the tolerance. Recall that we generate 2 10 6 events. 



This corresponds to the asymmetry A 1 = rjx/2 = 0.111 and the amplitude of cos0 to 
be ttA 1 /2 = 0.175. This is confirmed by the fit in frame M, see Table |3|. We would 
need 3 Nm rs 730 events to measure it with 3a significance in frame M. In frame F 
the asymmetry in this case is A 1 = —0.035 and requires Np ~ 7350 events for it to 
be measured with 3a significance. Thus, one needs at least max(iVM, Np) = 7350 
events to confirm the spin of t-quark to be at least | with unpolarised beams. In 
this case where the beams are not polarised, the asymmetries are smaller in frame 
F however the analysis in this frame does confirm that no new modulation has been 
missed, and thus reconfirms the spin-1/2 nature of the top. 

In order to improve the sensitivity, one might consider the case of polarised e + e~ 
beams to produce top quarks with larger polarisation. For example, 

(P e -,P e +) = (+0.80,-0.60) : 771 = —0.505, r] 2 = 0.000, rj 3 = +0.554 , 

which corresponds to much larger polarisation and hence a larger asymmetry A 1 = 
—0.253 in frame M. This requires only Nm ~ 140 events to measure A 1 with 3a 
significance. The azimuthal distribution for this beam polarisation is shown in Fig. ^ 
in both the frames M and F. In frame F however the asymmetry A 1 is smaller, 
see Table |3], hence we need a larger number of events, Np w 410, to measure it 
with 3a significance. Thus, one needs max^^f, Np) = 410 events to confirm the 
spin of t-quark to be at least | with this choice of beam polarisation, which is a 
large improvement over the unpolarised case. To rule out higher asymmetries with 
a higher degree of significance one still needs a larger number of events than this. 

3 Number of events required: N — f 2 /(A^) 2 , where / is the degree of statistical significance. 
Numbers with / = 3, for 3ct significance, are quoted. 
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Figure 4: The azimuthal distribution of lepton from decay of t-quark is plotted for 
(P e -,P e +) = (+0.80, —0.60) in frame M (top) and in frame F (below) using 2 x 10 6 events 
at partonic level (histogram). The best fit (green/grey line) to F^{(p) leads to the coefficient 
of the cos 4> modulation to be non-zero (given in Table ||) and all other modulations are 
absent in both the frames indicating the spin of t-quark to be ^. 



Next we discuss the case when the transverse polarisation of t-quark, 771, is zero. 
We arrange this by tuning the beam polarisations to appropriate values. This leads 
to A 1 pa 0, and hence in frame M the fit gives ai/ao ~ 0, and a flat distribution as 
shown in Fig. [5|. The top polarisations in this case are given as 

(P e - ,P e +) = (+0.792, +0.60) : ^ = 0.000, r} 2 = 0.000, 773 = +0.080 . 

We note that the longitudinal polarisation of t-quark, 7/3, though small is not zero 
and hence in frame F this leads to a non-zero value of the asymmetry A 1 and the 
cos0 modulation as seen in Fig. |^. In this case Np ~ 2 10 4 events are required to 
measure this asymmetry at 3a significance. This example re-imposes the need for a 
second frame F in association with the helicity frame M to measure and re-confirm 
the spin of a particle. 

4.2 Spin— 1 case: W— boson 

The W boson analysis is a much better advocate for the need of frame F, beside 
frame M. In the process under consideration, the W-bosons are produced (almost) 
on-mass-shell as a decay product of t-quark. Since the coupling of PF-boson is chiral, 
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Figure 5: The azimuthal distribution same as in Fig. |] for (P e -,P e +) = (+0.792, +0.60). 
This is special case when the distribution is flat in frame M and frame F is needed for spin 
determination. 



they are produced with high polarisation even in the decay of unpolarised top quarks. 
For the same set of events as used for the case of top quarks, the vector polarisations 
of the VF-boson are given by 



(P e -,P e +) = (+0.80, -0.60) : p x = +0.355, ^ = 0.000, 



Pz 



0.000 



From Eq.( |3.4| ) we know that the coefficient of the cos0 modulation is proportional 
to p x and non-zero in this case. The coefficient of cos 2<p modulation is proportional 
to tensor polarisation (T xx — T yy ), which happens to be zero 4 for this process in the 
helicity frame M. This leads to only cos 4> modulation of the azimuthal distribution 
as seen in Fig. ^| for frame M. Using the helicity amplitudes given in Appendix |C.1| , 
one can write the production density matrix for VF-boson, which is produced in 
the decay of t-quark, and we easily see that piy(+i, — 1) = pw(— 1? +1) = in the 
helicity frame M due to angular momentum conservation. Here, a higher spin (spin- 
1) particle is produced in the decay process of lower spin (spin-|), thus it can not 
span all its helicity states for fixed helicities of other particles and hence leads to 
Pvy(±l, +1) = 0. This is proven for the general case in Appendix |B] in the helicity 



4 We note that the asymmetry A 2 is zero in the helicity frame M for on-shell W bosons, but 
numerically we find it to be small but non-zero as the decay width of W is not very small and there 
is a non-negligible contribution from off-shell Ws. 
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Figure 6: The azimuthal distribution of lepton from decay of VF-boson is plotted for 
(P e -,P e +) = (+0.80, —0.60) in frame M (top) and in frame F (below) using 2 x 10 6 events 
at partonic level (histogram). The best fit (green/grey line) to F±{4>) leads to the coefficient 
of the cos (ft and cos 2</3 modulation to be non-zero in frame F (given in Table |j) and all 
other modulations are absent in both the frames indicating the spin of VF-boson to be 1. 
The red (dark grey) line show only the cos <j> modulation of the distribution. 



frame. However, in the boosted frame F the asymmetry A 2 measures T^ x — T^ y , 



which is non-zero in general, see Eq.( |3.14j ). In frame F we find A 2 = —0.054 which 



leads to a cos 20 modulation of the azimuthal angle in this frame, see Fig. ^. Here 
Np = 3100 events will be required to measure A 2 with 3a significance. Further, all 
the higher A J s (A 3>2 ) are found to be zero in both frames proving that the particle 
under consideration to be spin-1 and its production process to be CP- conserving. 
The asymmetries and fit parameters are listed in Table (| for both the frames. 

Next we look at a case where the azimuthal distribution in the helicity frame M 
is flat which would wrongly suggest that the particle is a scalar. The various vector 
polarisations are given as 

(P e -,P e +) = (+0.75, +0.60) : p x = +0.000, p y = 0.000, p z = 0.193 . 

In the helicity frame M, the asymmetry A 2 is zero due to the angular momentum 
conservation and A 1 is zero because it is proportional to p x , which is zero for the 
chosen initial state beam polarisations. The H^-boson appears to be spin-0 in this 
frame M with this particular beam polarisations. The asymmetries A 1 and the fit 



26 



(P e -,P e+ ) 


Quantities 


Frame M 


Frame F 


Reference 




A 1 


-0.266 


-0.038 






A 2 


w 


-0.054 


Fig.| 


(+0.80,-0.60) 


d/ao 






-0.418 


-0.059 




a 2 /a 


« 


-0.086 






A 1 


« 


-0.093 






A 2 


« 


-0.026 


Fig-0 


(+0.75, +0.60) 


ax /a 






« 


-0.147 




a 2 /a 


« 


-0.041 





Table 4: The table of asymmetries A 1 & ^4 2 and the fit parameter ai/ao & 02/00 for 
lepton's distribution from VF-boson decay for different initial state polarisations in frames 
M and F. We have the relation ai/ao = irA i /2, which is also observed numerically within 
tolerance (+10 -3 ). 



parameters aj/a are listed in Table £| for this case and the corresponding azimuthal 
distributions are plotted in Fig. [7[ Changing over to frame F leads to non-zero values 
of both A 1 and A 2 , see Table J|, and the corresponding azimuthal distribution visibly 
has the cos 20 modulation, Fig. |7|. In this case Np « 1.3 x 10 4 events are required to 
measure A 2 with 3a significance. This is the best example of a case where one needs 
a frame other than the helicity frame to confirm the spin of the particle, which is 
polarised with p z 7^ and T zz 7^ 0. We, however, note that this process is not the 
best process to study the spin of W-boson. For this purpose one should look at the 



pair production process e + e — > W + W as discussed in Ref. |p4 



Thus we conclude that one needs two different reference frames to measure and 
re-confirm the spin of a particle using the same set of events. Further, if the event 
set includes a cascade decay, one can construct the asymmetries for different parti- 
cles using the spin vector Sj for different particles and using the momentum ps of 
different final state particles. For example, in the above case, we could have used the 
momentum of the 6-quark to construct the asymmetries in place of the leptons. In 
the events with hadronic decay of Ws, one could use either of the jets to construct 
the correlators and hence the asymmetries. Thus, using different final state particles, 
we can find a larger set of asymmetries to re-confirm the spin of a particle, however 
we can not improve the significance of the measurement by combining different cor- 
relators for the same set of events. A larger event sample is necessary to improve the 
statistical significance of the measurements. 
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Figure 7: The azimuthal distribution same as in Fig. |6| for (P e -,P e +) = (+0.75, +0.60). 
This is special case where frame M has flat distribution and frame F is needed for the spin 
determination. 



5. Discussions and Conclusions 

In this paper, we constructed observables to measure the spin of a heavy unstable 
particle produced at a collider by harvesting the spin dependence of the azimuthal 
distribution through the quantum interference between the different helicity states, 
i. e the non diagonal elements of the helicity density matrix. The aim is to construct 
observables that are sensitive to the highest rank-2s tensor polarisation of a particle 
with spin-s that lead to a cos 2s(p modulation of the azimuthal distribution. Such a 
method has been known for a long time and we have provided an analytical under- 
standing of it. In particular, the novelty of our approach is the construction of two 
reference frames where in one of them the spin basis is subjected to a Wick rotation. 
The latter mixes the longitudinal polarisation and the transverse polarisation in the 
production plane and therefore the spin modulation in the azimuthal angle is sensi- 
tive to this mixture whereas in the standard approach the longitudinal polarisation 
is integrated away and does not contribute to the usual azimuthal asymmetries. The 
construction of two frames allows within the same experiment and with the same 
event sample to cross check the spin measurement based on azimuthal asymmetries. 
In some cases this can be crucial since the usual transverse polarisation tensor/vector 
can be zero either accidentally or for a dynamical reason and therefore would lead 
to a wrong conclusion. This can be rescued in the second frame provided the longi- 
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tudinal polarisation is not zero as well. We have shown examples in the decay of the 
top and the W where this occurs with SM production and decay mechanisms. In the 
appendices we consider more general couplings and decays than those that describe 
the SM particles. This helps in obtaining a set of conditions on the production and 
decay mechanisms for some of the asymmetries to be non-zero and hence the spin 
to be measured. One drawback of the method however, as outlined in the present 
paper, is that it requires complete reconstruction of the test particle's momentum 
which is necessary to build up the needed spin vectors. If there are too many invisible 
particles this might not be possible especially in a machine like the LHC where the 
partonic centre of mass energy is not fixed. The same drawback also affects other 
methods of spin reconstruction. We feel however that is worth investigating how the 
method we have described can be exploited in combination with other methods or 
by making some mild assumption on the spectrum of the event or the underlying 
physics. 
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A. Rotation matrices d J mn (0) 



The general form of the d function is given in Eq.( |1.2j ). For completeness and al- 
though these can be easily found in many textbooks, we explicitly write the d func- 
tion up to spin-2. To avoid clutter we take as short-hand notation c = cos(#/2) and 
s = sin(0/2) then all the d{ nn useful to our study are given below. 

1 



J 



: d° 



I ■ J? 





c 2 


-V2cs 


s 2 




V2cs 


2c 2 -1 


-V2cs 




s 2 


V2cs 


c 2 
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3 

d 2 



Vise 2 
VSs 2 c 



1 = 2 

c 4 
2sc 3 

2s 3 c 



d 2 



-V3sc 2 
c - 3s 2 c 
-s + 3sc 2 
V3s 2 c 



-2sc 3 
-3c 2 + 4c 4 
■a/6sc(s 2 - c 2 ) 
3s 2 - 4s 4 
2s 3 c 



y/Ss 2 c 
3sc 2 
3s 2 c 
V3sc 2 



s 
c 



— S J 

V3s 2 c 
-V3sc 2 



\/Qsc{s 2 - c 2 



1 - 6c 2 + 6c 4 
— t/Gsc(s 2 — c 2 ) 
VQs 2 c 2 



-2s 3 c 
3s 2 - 4s 4 
/Qsc(s 2 - (?) 
-3c 2 + 4c 4 
2sc 3 



s 

-2s 3 c 
3s 2 c 2 
-2sc 3 



VQs 2 - 2 



B. Decay density matrix for higher spin particle 

As a short-hand notation we now define C = cos 6 and S = sin9 which enter the 
expressions for the density matrices of higher spin particles, namely s = 3/2 and s = 2 
briefly discussed in the main text. The corresponding normalised decay matrices are 
calculated from Eq. ^/F] and using the explicit expressions for the d matrices. 



B.l Spin-| particle 

For the decay ||, I) — > \si, h) + |s 2 , h), the decay density matrix is given by 



ri 


, 3 

(+ 2' 


+ 2> 


_ (1 + 2 7l ) + 3(«! + a 2 )C + 3(1 - 2 7l )C 2 + (a 2 - Sc^C 3 

8 


ri 


, 3 


+ 2> 


_V3S[ ( ai + a 2 ) + 2(1 - 2 7l )C + (a 2 - 3 ai )C 2 ] w 

8 


r§ 


r 3 

:+ 2' 


1 

T 


_ v/3 S 2 [ (1 - 2 7 i) + (a 2 - 3ai)C7 ] w 
8 


ri 


, 3 

;+ 2' 


3 

~r 


(«2 - 3ai)^ 3 i30 
8 


ri 


, 1 

(+ 2' 


3 X 

+ 2> 


_ 73 5 [ ( ai + a 2 ) + 2(1 - 2 7l )C + (ot 2 - 3 ai )C 2 ] ^ 

8 


ri 


, 1 

(+ 2' 


+ 2> 


(3 - 2 7l ) + (3a 2 - 5ai)C - 3(1 - 2 7l )C 2 - 3(a 2 - 3ai)C 3 ] 

8 


ri 


, 1 

;+ 2' 


1 

~r 


S [ (3a 2 - on) ~ 3 ("2 - 3a!)C 2 ] i(j> 
8 6 


ri 


, 1 

;+ 2' 


3 


V3S 2 [(1- 2 7l ) - (a 2 - 3ai)C7 ] m 
8 


ri 


, 1 

~2' 


+ 2 } 


V3 5 2 [(l- 2 7l ) + (a 2 - 3ai)C7 ] ^ 
8 6 
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r§ 




r§ 


' 2' 2 J 


r§ 


* 2' 2 J 


ri 


, o o . 


r§ 


,3 1. 
-2' + 2 } 


ri 


_3 _1 

' 2' 2 J 


ri 


, 3 3 
~2'T 



where, 



and 



S [ (3a 2 - ai) - 3(a 2 - 3cti)C 2 ] _ i0 

8 6 
(3 - 2 7 i) - (3« 2 - 5ai)C - 3(1 - 2 7 i)C 2 + 3(a 2 - 3«i)C 3 

8 

>/3 5 [ («i + a 2 ) - 2(1 - 2 7l )C + (a 2 - 3 ai )C 2 } J<p 



e 



(g 2 ~ 3ai)^ 3 _ i3( ^ 
8 

VS S 2 [ (1 - 2 7 i) - (a 2 - 3ai)C ] „_ i20 



e 



3 5 [ (ai + a 2 ) - 2(1 - 2 7l )C + (a 2 - 3ai)C 2 ] 



e 



(1 + 2 7l ) - 3(ai + a 2 )C + 3(1 - 2 7l )C 2 - (a 2 - 3a 1 )C 3 



3/2 
°3/2 ~~ 




\h\ 


< Si, 




< s 2 


3/2 
fl l/2 — 




\h\ 


< Si, 




< S2 


3/2 
fl -l/2 ~~ 




\h\ 


< Si, 


H. + 5I 


< s 2 


3/2 
°-3/2 — 


h 


\h\ 


< Si, 


It + fl 


< S2 



B.2 Spin-2 particle 

For the decay |2, /) — * \si, h) + |s 2 , / 2 ), the decay density matrix is given by 

T 2 (+2, +2) =[A + AA 1 C + QA 2 C 2 + AA 3 C 3 + A 4 C A ] 
T 2 (+2, +1) = 2 [ A 1 + 3A 2 C + 3A 3 C 2 + A 4 C 3 ] S 
T 2 (+2, +0) = V6 [ A 2 + 2A 3 C + A 4 C 2 ] S 2 ^ 
T 2 (+2, -1) = 2 [ A 3 + A 4 C ] S 3 
T 2 (+2,-2) = A A S 4 e 1 ^ 



(B.l) 



3/2 _ 3/2 3/2 _ 3/2 3/2 3/2 

EzV EjOj Eiflj 



(B.3) 
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T 2 (+l, +2) = 2 [ Ai + 3A 2 C + 3A 3 C 2 + A 4 C 3 ] S e _< * 

T 2 (+l, +1) = 4 [ 1 + 2(Ai - 3/3) C - 3A 2 C 2 - 2A 3 C 3 - A 4 C 4 ] 

T 2 (+l, +0) = 2V6 [ 2(/3 + A 2 C) + S 2 (A 3 + A 4 C) ] S e i(t> 

T 2 (+l, -1) = 4 [ 3A 2 + A 4 S 2 ] S 2 e i2 * 

T 2 (+l, -2) = 2 [ A 3 - A 4 C ] S 3 

T 2 (+0, +2) = V6 [ A 2 + 2A 3 C + A 4 C 2 ] S 2 e~ i2<t> 

T 2 (+0, +1) = 2VQ [ 2(/3 + A 2 C) + S 2 (A 3 + A A C) ] S e"** 

T 2 (+0, +0) = 4 [ 45 + 3A 2 S 2 + 3A 4 S 4 ] 

T 2 (+0, -1) = 2VQ [ 2(/3 - A 2 C) + S 2 (A 3 - A 4 C) ] S 

T 2 (+0, -2) = V6[A 2 - 2A 3 C + A 4 C 2 ] S 2 e <2 * 

r 2 (-l, +2) = 2 [ A 3 + A 4 C ] S 3 e" 13 * 

T 2 (-l, +1) = 4 [ 3A 2 + A 4 S 2 ] S 2 e- l2<t> 

r 2 (-l, +0) = 2V6 [ 2(/3 - A 2 C) + S 2 (A 3 - A A C) ] S e-** 

T 2 (-l, -1) = 4 [ 1 - 2{A 1 - 3/3) C - 3A 2 C 2 + 2A 3 C 3 - A 4 C 4 ] 

T 2 (-l, -2) = 2 [ A l - 3A 2 C + 3A 3 C 2 - A 4 C 3 ] S e* 

T 2 (-2,+2) = A 4 S 4 e-* 4 * 

T 2 (-2, +1) = 2 [ A 3 - A 4 C } S 3 e'^ 

T 2 (-2, +0) = VQ [ A 2 - 2A 3 C + A 4 C 2 ] S 2 e _i2 * 

T 2 (-2, -l)=2[A 1 - 3A 2 C + 3A 3 C 2 - A 4 C 3 ] S e~* 

T 2 (-2, -2) = [Aq- AA X C + 6A 2 C 2 - 4A 3 C 3 + A 4 C 4 ] (B.4) 

where, 



A = 


a 2 + 4a 2 + Qal + 4a 2 x + a 2 2 


16 Ez 




A l = 


a| + 2a 2 - 2a 2 x 


-a 2 2 


16 E* a? 




A 2 = 


a 2 — 2a 2 , + a 2 2 




16 $>? ' 




A 3 = 


a 2 . - 2a\ + 2a 2 x 


-a 2 2 


16 J2i a 2 




A 4 = 


a\ — 4a 2 + 6a 2 , - 


- 4a 2 1 + a 2 2 




a 2 


P = 


aj - a 2 ! 


a 2 


EX ' ' = 


EX 



(B.5) 



^ = iZ)l- M Wi-2l 2 I'll < |ii-2|<a 2 
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a l 



Cln = 



a_, = 



a-2 



Aix 

h 



TT^I-^LJ 2 |Zi| <minsi,s 2 



47T 

il 



47T 



° Y.\ M l,h+2? \k\<si, \h + 2\<s 2 . (B.6) 



47T 

h 



C. Helicity amplitudes and the analysing power 

In this section we give expressions for the helicity amplitudes pertaining to 2-body 
decay processes of spin-| and spin-1 particles. The expressions are derived for a 
general dimension-4 effective operator describing the coupling of the particles. This 
will permit to give the different analysing power coefficients. 

We take the mass of the mother particle to be m and that of daughters to be m\ 
and m 2 , the polar and azimuthal angle belongs to the first particle with mass m±. 
The energy and the momentum of the daughter particles are given as 

m 2 + m\ — m\ m 2 + m\ — m\ 



Ei = , E< 



2m 2m 
_ a/ ((m + m 2 ) 2 - m\)({m + mi) 2 - rnf) , . 

from 2-body decay kinematics. Below we discuss the 2-body decay of a fermion and 
a vector boson into two massive particles. 

C.l Decay: ||, A) — > ||, Ai) + j 1, A 2 > 

For this decay the helicity for the fermion A = ±1/2 will be denoted as A = ±1/2 
and for the bosons A = ±1 as A = ± such that the helicity M(A, Ai,A 2 ) writes as 
M(±,±,±) = M(±|,±|,±1). 

The decay vertex it taken to be /i7 M (ClPl + Cr-Pr),/^^ with real Cl,r and the 
amplitudes are listed below in the rest frame of the decaying particle: 

M(±,±,±)= [-(C L Pr - C R P?)] (-sin?) 

M(+, +, 0) = [-(C L PfP 2 - - C R P+P+)] (+ cos 

M(+,+,-) = 
M(±,-,+) = 

M(±, 0) = [+(C L P+P 2 + - C R PfP 2 -)] (- sin Q 
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M(±,-,-) = [+(C L P+-C R P 1 -)} (+cos0 



cos • 



M(- +,+)= [-(C7 L Pf - C7 fl P+)] e— ( + 

M(-,+,0) = [-(C L PfP 2 - - C7 K P+P+)] (+sm£) 

M(- +,-) = 
M(-,-,±) = 

M(- -,0) = [ + (C7 L P+P+-C7 K P-P-)] ( 
M(-, -, -) = [+(C7 L P+ - Pf)] f + 



cos 



sin 



(C.2) 



Here, the terms in the square brackets are the reduced matrix elements, (M. s Xlj \ 2 / V^rr) 
and the d s X Xl _ X2 functions are enclosed in round brackets. The symbols P^ 2 are de- 
fined as 

Ei + mi ± p 



Pi 



± 



m 



y/E{ + mi 



,± = J_ E 2 ±p 
2 \/2V E 2 TP 



(C.3) 



Using the expressions of the reduced matrix elements, the analysing power a for this 
decay can be written as 



a 



(C% - Cp{\ -xj- 2sjVj + (g - gg - 2(xg +g 
(C| + C 2 L )(1 - 2xj + x 2 2 + x\x\ +x\- 2x1) - \2C L C R xix\ 



2 ' 



(C.4) 



where Xi = rrii/m. For the decay of top quark, t — * bW, with m\ = = within 
the SM we have C R = leading to a = -(1 - 2x|)/(l + 2x|) 0.38. 

C.2 Decay: |±,A) -> |±,Ai) + |0,0) 

As done in the previous section the helicity for the fermions A = ±1/2 will be denoted 
as A = ±1/2 such that the helicity M(A, Ai writes as M(±, +) for M(±|, + 

For this decay the helicity amplitudes, M(A, Ai) = M(±, +) = M{+\, +\). The 
decay vertex it taken to be (C L P L + C R P R )f 2 S with complex C L>R and all the 
amplitudes are listed below: 



M(±,±) = 
M(±,-) = 
M(-,±) = 







■ + c L 


Pf 






v/2 






J? 


"±C L 


Pi 






y/2 






A~ 


~ + c L 


Pf 



Vi<t> 



e 2 



+»</ 
e 2 





±cos- 



-sm- 



V2 



± sin 
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M(- 



c R p+ + c L pr 

V2 



— i<jj 



+ cos ■ 



(C.5) 



Here, the reduced matrix elements, (Ai s x /v2tt), are given in square brackets and 
the d s XXl functions in the round brackets. The symbols Pf" are same as in Eq. (|C.3|) . 
Using the expressions of reduced matrix elements we get the expression for a, the 
analysing power of the spin-| particle, as 



a 



~(\C R 


2 _ 


|C L | 2 )v/l + (x?-x|)2-2(x? + x|) 


(\Cr\ 


2 +l 


Cl\ 





(C.6) 



where = rrii/m. Thus we need the scalar to have parity violating couplings, 
\Cl\ 7^ \Cr\, for the analysing power to be non-zero. However, for a neutral scalar, say 
the neutral Higgs boson of the MSSM, we have Cl = C* R in other words \Cj\ = \Cr\ 
leading to a = 0. The same occurs in a CP conserving MSSM with any of the neutral 
Higgs boson. Thus we should chose processes involving squarks for spin measurement 
in the decay of gauginos. 

C.3 Decay: |1,A) -> |§, A x ) + |§, A 2 ) 

We take the same convention as in |C.1| with the same operator for the interaction. 
We find 



M(+ 
M(+ 
M(+ 
M(+ 
M(0 
M(0 
M(0 
M(0 
M(- 
M(- 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+i 



—i 



+i 



+i 



—i- 



+i 



■ C R pjp 2 + C L p 1 pj 

2 

Cr pXpX + C L plp^ 
V2 

Cr P1P2 + c l pXpX 
V2 

Cr P1P2 + C L p\p2 
2 

■ Crp\P2 +Cl Pi pX 

2 

Cr pXpX + C L P1P2 
V2 

Cr P1P2 + C L pXpX 
V2 

Cr P\p\ + C L pXp 2 
2 

■ CrpXp2 +c l p\p\ 

2 

Cr pXpX + C L p\p 2 
V2 



I COS I 



f —= sin 9 

\V2 



2d 
COS - 

sin - 



-1 

7! 



sin # 



— = sin 6 1 
Vv/2 

— = sin 
\/2 



cos( 



-i<f> 



+ 1 " fl 

— = sin U 

V2 



sm 
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M{-,-,+) 
M{-,-,-) 



—i 



+i 



Cr y x p 2 + C L p\p\ 
Cr PiPt + C L p\p 2 



cos 



-i<f> 



—= sin 9 

V2 



(C.7) 



In Eq. |C.7| the terms in the square bracket stand for the reduced matrix elements 
(M. Xl \ 2 \J^)i the terms in the round brackets are the d s XXi _ X2 functions. The sym- 
bols pf 2 are defined as 



Pi 



± _ E 1 + m 1 ±p 
y/Ei + mi 



P2 



E 2 + m 2 ±p 
y/E 2 + m 2 



(at 



Using the expressions for af and reduced matrix element we get expressions for two 
parameter a and 5 as 



a 



\2C L C R x x x 2 + (C| + C\) [2 - (x 2 - x|) 2 + {x{ + x 2 )] ' 

4C L C RXl x 2 + (Cl + Cl)[{x\ + xf) - (xj - x 2 ) 2 ] 
\2C L C R x x x 2 + (C 2 + C£) [2 - (xf - x 2 ) 2 + (x 2 + x 2 )] ' 



(C.9) 
(CIO) 



where Xj —rrii/m. 

If the final state fermions are massless, xi — > 0,X2 — > 0, one obtains a 
c2 l)/(C r + C£) and <5 — > 0. This is the case for the decay of W and Z bosons into 
massless fermions. Further, for the decay of Ws, within the SM we have Cr = 
hence a = — 1. 



C.4 Decay: |1,A) -> |l,Ai) + |0,0) 

For this decay the helicity amplitudes, M(A, Ai) = M(+, +) = M(+l,+l). The 
decay vertex it taken to be CvvsQ^uV^Vi with real Cyvs and the helicity amplitudes 
are given by: 



M(+,+) 
M(+,0) 

M(+,-) 
M(0,+) 
M(0,0) 
M(0,-) 



-Cvvs] e +J ' 



mi 



2^ 
COS - 



sin # 



v/2 



-CWs] e +l * { sin 2 r 



-C 



vvs\ 



sin 6 1 



—^vvs — 



mi 



-C 



vvsl 



(cos 6) 
— sin 8" 
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M(-,+) = [-Cws] e~» (sin 2 ^ 
M(-,0) = 



mi 



sin # 



v/2 



M(-, -) = [-C vvs \ e- u 
This leads to a\ = hence a = 0. 5 is given by 



cos 2 - 
2 



(C.ll) 



l + (x 2 -:r 2 ) 2 + 2(5x 2 -x 2 ) 

Such decays occur the models of extra-dimensions, for example, 
Z W etc. 



(C.12) 



C.5 Decay: |1,A) -> |l,Ai) + |1,A 2 ) 

For this decay the helicity amplitudes, M(A, Ai, A 2 ) = M(+, +, +) = M(+l, +1, +1). 
The decay vertex is taken to be CyyyT^pV^V^V^ where C V vv is real. We only 
assume here a standard gauge tri-linear coupling 

T„ U p = [9^u(q - Pi) P + 9v P {v\ - V2)p, + 9 P p(P2 - q)A 

with q being the 4-momentum of the mother particle and p\^ is the 4-momentum of 
the daughter particles. All momenta are assumed incoming at the interaction vertex. 
With this notation the non-zero helicity amplitudes are given by: 



M(+,0,0 

M(+,-, 
M(0,+,+ 
M(0,0,0 
M(0,-,- 

M(-,+,+ 

M(-,0,0 
M(- - - 



[2C vvvP ] e +irt > 
-2C VVV p 



- sin 9 

7T 

E X E 2 + p 2 



17111712 



— sin 9 

~7T 



[2C vvvP ] e +i<t > 
[2Cvvvp] (cosO) 



- sin 9 

7T 



-2C VVV p 



E 1 E 2 +p 2 



i7iim 2 



(cos 9) 



[2Cvvv p] (cos i 



\2Cvvvp] e 

-2C VVV p 

[2C vvv p] e 



sin 9 



V2 
E^ + p 2 ' 
mim 2 
— sin# 



-i(f> 



— sin 9 



V2 



(C.13) 
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This leads to a\ — a\_ x — and hence a = and 5 = 1. Example of such decays, in 
the models of extra-dimensions, are W {1] -> W {0) Z^°\ W {1) -> etc. 

Above we saw that the parameters a and 5 have a simple expressions in terms 
of masses and couplings of the particles involved. This can be simply added to a 
spectrum generation code, such as SOFTSUSY f3§], SuSpect [|36|, SPheno |7| for SUSY 
models, and one can quickly know which decay channel is the best for estimation of 
the particle's spin. 

D. Higher spin particle disguising as lower spin particle 

If a higher spin particle is produced as a decay product of the lower spin particle, its 
spin orientations are restricted. This makes the particle appear as of a lower spin in 
frame M see Sec. [0|. Since the total differential rate is the product of production 
and decay density matrices, 

dv = ^jP s (U') x V s (1,1') rf$ n . 

Thus, for the decay distribution to have a 2s0 modulation, we must have p s (s, —s) = 
p s *(—s,s) ytz o. Now if this spin s particle were produced in the decay reaction 
\j, m) — > \si, li) + \s,l), then the production density matrix is given by 

P s (M') = E M Z M U 

Thus we have extreme off-diagonal term given by 

P S ^ s) = (^P) E d Us d U + s MlMtr (D.2) 

^ ' m,h 

For this to be non-zero, we must have 

\h — s\ < j and |/i + s\ < j 

for at least one value of l\. However, this condition is never satisfied (for any l\) 
when we have s > j, i.e. p s (s, —s) = for s > j. This leads to the absence of the 
highest mode in the <fi distribution in frame M. This is numerically demonstrated 
for the sample of W boson production from the decay of t-quark. Since the helicities 
are invariant only under the boost along the momentum (which does not changes 
the direction of the momentum), the density matrix goes through a similarity trans- 
formation when boosted in any other direction. Thus in frame F, in general one can 
have a non-zero value for p s (s, —s) and hence the 2s0 modulation of the azimuthal 
distribution. 



-m 
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